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This investizati nun was conducted to deterrmine the effects of 
anguler base motion interfererce on a geometric stabilizeti n system 
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CHAPTER I 


IhTRODUCTION 


Recent scientific and e: ginesring advances have added satellites 
Shismeatitewic messes «oO lune: kinds: of, veniel>s used by man... ine 


successful us2 of these vehicles recuires 2 rears of cortrojJling their 
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propulsicn. The onerating region of these vehicles and the  ro:ulsior 
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vector with the instantaneous velocity of the vehicle. The task of the 
control system is to identify the vehicle velocity vector and maintain 
its magnitude and orientation in a prescribed manner. Most important, 
the control system must cause the vehicle to attain the predetermined 
value of the velocity vector at thrust cut-off. This means that the 
purpose of the control system is to control the magnitude and direction 
of the velocity of the vehicle at thrust cut-off. 

In order to control the orientation of the vehiclels PietOCley 
vector, the control system must (1) establish and maintain a reference 
orientation in some usable form (2) provide a means for comparing the 
actual velocity vector orientation with the reference orientation and 
(3) provide a means for altering the velocity vector to make the dev- 
jation from the reference become zero. The reference orientation need 
not be constant during the flight. It may coincide with the desired 
flight path or with any other prescribed path. It is possible to dis- 
tinguish btetween an inertial and an instantaneous reference orientation. 
Themconcro:-systen's function: is to bring the velocity vector into 
coincidence with the instantaneous reference orientation. 

It is desireable from the military standpoint to eliminate any 
dependence of the vehicle on external sources of information, There 
are systems that are not self-contained, requiring some of the components 
of the control system to be located on the ground or elsewhere outside 
of the vehicle. There are other systems that require external radiation 
or external reference points to establish the reference orientation. 
Inertial guidance systems have been constructed that are completely 


independent of external sources of focnacsomt These systems use 


* Superscript numerals refer to similarly numbered references in the 
Bibliography at the end of this paper. 


(3) 





measurements made with respect to inertial space to establish the 
reference orientation and to measure deviations from the reference, 

Since the control system depends on a process of relating the v2locil, 
vector of the vehicle to the reference orientation (now assumed to be in 
inertial space) any uncertainty in the repr2sentation of the reference, 
or in the information that corresponds to the representation, is equivalen. 
to an wicertainty in the orientation of the velocity vector. The ability 
to control the orientation of the vehicle's velocit: is impaired to the 
extent that the reference orientation is uncertain. To obtain high 
accuracy in controlling the vehicie trejector;, it is therefore necessary 
to maintain the reference orientation witl hich accuracy. 

The most elementary task of the control syst-m is to maintain the 
reference orientation when this orientation is fixed in space. Ge2orm=stric 


4 eye e he £ = 
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stabilization (2) is the ability of the sjsten 

orientaticr: in the presence of 2xterral sources cf interference, ‘The 

ability to track a changing reference orientetion and to command the veloc: 
‘ 


weCLOr LO aecuire 4 prescribed “OrTlencztion with rescect to the instant— 


a 


anecus reference is the guidance function of the cortrol system. However 
a measure of the capability of the control system is 1ts succ2ss in 
achieving geumetrical stabilization. 

Such a stabilizstion system has be2n instrumented using three 
Single degree of fre2dom integrating gyros mounted cn a gimbaled platforn 
with their input axes along thr22 mutually verrendicular directicns. A 
line schematic diagram of a three axis zimbal system used for supporting 
sucn a platform is shown in Fig. eed The output sicnal of each gyro 
is proportional to the integral of the angular velocity of the platform 
about the respective axis. The outputs of these gyros actuate 


gimbal servo drive motors. These cause the platform to rotate so as to 
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a. THE ELECTRICAL POWER SUPPLIES, ELECTRONIC UNITS, COMPUTERS, '‘CCNNECTIONS, 
RACKS AND OTHER COMPONENTS NECESSARY TO CCMPLETE AN INERTIAL GUIDANCE 
SY OLEM AME NOT REPRESENT EDIIN 7 HIS. FIGURE. 


3. THIS ILLUSTRATION IS BASED ON FIG. 4 OF WRIGLEY, WOODBURY AND HOVCRKA (2) 
AND FIG. 6 OF DRAPER AND WOODBURY (4). 


186 1-1 Line schematic diagram showing essential mechanical elerents of inertial guidance system 
based on rotation of the inertial reference package with the indicated vertical. 





null the gyro signals. At null indication of the zrro the platform is 
essentially restored to its original orientation in inertial space, 
Figure 1-2 is a pictorial diagram of one such instrumentation. Such a 
system may be considered as three, single axis space integrators (2) 
operating in parallel. 

The essential features of a single axis svac2 intevrator are shown 
in the block diagram of Figure 1-3. This investisation was conc2med 
only with the geometric stabilization problem. For the remainder of 
this paper it will be assured that the command signal is z2ro, and the 
actuating input to the system is an interfering torcue acting on the 
controlled memcrer. The performence of a systen. }.ossessing g2ometric 
stabilization may b2 lescribed by the following steps: (1) The inter- 
fering torque acting on the controlled member causes an engular rotation 
with respect to inertial space about the gyro input axis. (2) 7 

ngular rotation caus2s a orecessiorn of the gyro with r2spect to the 

controlled member, producir.g a sigmel from the sienal generator, which 
is proportional to the anzle of precession. (3) The gyro signal is fe? to 
the servo motor which them applies a torque to rotate the contrelied 
member backvto its initial orientation, nulline the cyro oul: ut si ma, 
System demping is provided by comrensation netwerks., 

The basic single axis system than, cc:.sists of 4 2zyro nackess 
mounted on a controlled member and 2 torcu2? gancrstor. The r:ference 
coordinate system is fixed in the cor.trclled membor 
the controlled member orientatior. mainteins the reference coordinate frame 

With the reference coordinate systen estatlishei within the nissile 
the means of controliinyg the missile to its desired ortertation within 
that frame may b2 cor sidered. It is through control over the orient- 


ation of the missile, that control is ex-rcised over the vehicle's 
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This Diagram is »besea on Sig. 5 Of The Sherman !*%. YSairchild 
Publication “und Paper No. FF-13, Institute of the Aeronautical 
Sciences, New Yorx, Janusry 1955. (Used Vith Permission) 


Pictorial Diarran Of a Three-Axis Inertial Space Rotation Intesrator 
Basei On Three Single-Axis Inertial Space Rotation And Command Sig- 
nal Receivers Vitn Servo Drives For The Three Gimbal Axes. 
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velocity vector direction, In order to fur.ctiqn, such a control system 
must include a means of generating the desired orientatior of the vehicle 
in the sefersnce freme, a means cf measuring :.ctual missile orisntation, 
ani a means of ch:anginy it lo minimize or 2liminate any deviations 
between the dzsircd and actual orientatio:,. Oince tiis investivatior. is 
concerned with stabilization, it will :2 essured thet the missile has 

the aviliey ol att-ivins the desired orier.t-ticm.ar.t the cunction 207 
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This Figure is based on Fig. 1-7 of "Inertial Juidance” (S)} a 
monograrh by C.Se Draper, Wewe wrigley, and Sidney yeas, mb- 
lished by tke Instrumentation Laboratory, “.I.7., Cambridge, 
Yassachusetts, Aucust 1957. (Used with permission.) 
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reference system uses servo motors as torque renzrators snd the missile 
control system uses gimbaled thrust engines for thet purvose. tlimiration 
of the functional duplication of the combined systems can result in 
saving weight, and possible simplification of the over-all] system. 
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up of elements involving the direction cosines reletins the body 
and inertial axes set. In general these elements will dzpend on 
the order in which rotations about the three inertial axis occur. 
However, it is shown in Appendix A that for rotaticns which are 


small in the sens2 that the product cf two angies is much less 
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than one radian, he expression (2-1) becomes, in matrix fora: 
(2-2) 


‘ 


I! oh A : ; 
(I-8). 1 (1-2 7 (25): | (1-2). 
a I ay 3 

| 

| 


” 
- 
! 
1 
ad 
we 
ra 
i 
-~ 
4 
/ 
So) 
bl 
(XN 
‘al 
( 
ee) 
~~ 
>< 


Mao ree cn, = Noscriacs > Ne eC rt ee ly ot Ce eae 1S. ore 


eneularewelocity A) and the integrals of thess commer us. (A) 


Lg M4 Ld - a tr -~ oo mw a _ ~ . : 
Mir me yl al ee odyiakos 196 hot. Veivy p be shit -aSicw 4 =) ie 
«a + ee + -~ + 
OW elms (C7201 vile Orie “cf “rotaticn. 
e ae e Nas rad re Ke ase 4 ae 
The expansion of ([=2) Sives for the cadror 1.45 of nage 


WeLOCtuy 2.07.7 Une inertia) exis: 


UnCerratitie tice eesrecsioan sives ths apts Veen ric phe 
bod yanas ur asdsobolt. the A<ineriie ) a%es: 


(2-4) 


—_— 
tea 
' 
C 
we 
~~ 
fi. 
ct 
l 
o~ 
+4 
t 
td 
ad 
L 
c+ 
o~ 
“af 
t 
bw 
— 
o-~ 
r-4 
' 
Se) 


7) 
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. = —~ “ va SG . _— + Pa 4 ~~ 
Sacew thom aso Sleeve ts “hich SE 2Oo ne orrOr S41 7ne) mrs seer tot 


sag a : ee foes ao re Se ae 2 Foe | 4 zie 
ees CeO EO Se CO Ay 6S se ee. ee Cae yp Og. Eo ee ees as 


fe) 
ce 
cfr 
yg 
() 


d 
ra 
pte 
To a 
J 
* 
D 
WW 
a8 
U 
i) 
aus 


CQ = -~ ” ~~ » 7 . + ~ 4 
emer ctr ene OCCU eh (shee fro e. Ae on 


omy . + = are 1s ; : = ees a ae See | $ 
Plat reC. 2o2C7S ees ALON 1S ees 1 ce ic ie 
- 4) Os 4° se $3 LSet sicae eae Re 4 ~ 
eo ee re on re 1 Oe OC oe re eo eee re ee a 
Vee ~ 4 i, x ales ea " ade 4 
Cee Oo eo oat) Ue eae tele ty CRIS 1 one ee Yaseen 


2 a e * é 
Oooo Cleat) ne Derlersanc> ef oro test 


(5) 


Robinson, 


The 


aD 


ouation (2-4) may be inter, reted from a sing acme Utes Wels? 


4U 


Mev Ont wiicr 2G Une  irtie Or Itscaling the netursict the errer 
in orientation which may develop. Surctose that there ar: tvo obser- 


Sion. ONC 2 ale ttiertiy ireterense rage Arder]= gn the Seca re sreace 


3 
«D 
<S 


On. Cr His j0rn Aas 


4) 


pe 


frame. Jach observer ccencent laces his attent 
ana wis sable to. recerd retetion of the beoiy with resrect te inertial 
Space about that axis. In addition the observer in body axes can 


Gavsc ere cauien Of whe body “th rescest: to inert-a) stace about ois 


Was eS ON et (UNS oct eS SUS jected to ar larcitrery rovation: the 


(18) 





obs@rver in inertial axes sxes 1 rotation about his X-axis which ts 

the sum of the terms o: the richt in (2-4) while the observer in 

body ax2s sees a rotation which is 2 ual to the first term on the 

risht in (2-4). Let the body now be brought to rest at some time (t,). 
Suppose now that the observer in hotly axes caus2s the body to rotate 
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for rotating the body are exactly this control system. The time 
period (0-t,) may be thought of as theo dynamic time lag in the system 
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error sipnals in the sy3uem.. Thisverror stynal itself recresents en 
error in orientation of the Mody with resrect to inertial space, but 

it is tr-nsient in nature. The system is aware of it and eventually, 
it will be driven out. However, there is in addition to this transtent 
error, the deviation expressed in (2-5) which the system is not aware 
of and hence, can do nothing about. The advantage of the expression 
(2-5) is that it express?s the error in orientation due to the effects 
of finite rotations ir !ependently of any transient d2viations in 
orientation due to system dynamics which may exist. 

Tne determination of the magnitude of the error in ori2ntation 
discussed above, and expressed in equation (2-5) is, at least in theory, 
straishtforwerd. The esuations (2-2) taken toyether with the ver- 
Pommnce ecudulens Ose <Sis.en tor set of equaticns Which specu, 
the motion of the missile 


Sowith, sosrect VCOrinervial Specs a. bot 


the inertial am body coordinate frames. It is only n2erssary to 


solve these ecuctions for the appropriate rotations and perform the 


subtraction indicated ty (2-5), and the error in oriznt-tion is 


we 


ebtained. PUntoOrcunately. Unsse Situations re mon=linear, so-vuh=t- ne 
solution in closed form for the sen2ral case is vossible. However, 
for a specific syst2m, solution by numerical methods is possible. 
Goerverore, if a System which is-more or less tyaicalef these likely 
to be ercountsred invpractice 15 specified. then useful in-ormation 
Concernias the magiiv of error likely to arise in mractical 
Syovemmunaer tle 17 ert ering ei fects of missilevanecular motion. micnv 
be obtained. The specification of such a system and the analysis 
indicated are the purpose of the remaind2r of this paper. 


Integratine Gyro Errors 


Before proc2ding with the specification of a system for analysis 
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one other error in orientation which can develcp in a system such as 
that shown in Fig. (2-1) should be discussed. This is essentially an 
instrument errcr occurring in the single degree of freedom integrating 
gyro. It is of interest because its effect on the system perfcrmance 
is, in some respects, remerkably lite the effect due to finite rotat- 
ions discussed above, This analysis will attem t to distinguisl. between 
the two effects and effect scme comparison of their relative magnitudes. 
For purposes of discussion, the geometric effect already described 


(5) 


while this new 


(7) 


error producing effect will be referred to as kinematic drift. 


will be referred to as the finite rotaticn effect 


The origin of kinematic drift lies in the manner in which a 
Single degree of freedom gyro performs its function. The principles 
of operation of these devices have been covered in detail Slesinere, ate 
so that a simole functional descrirtion of their operaticn will suffice. 
Fig. (2-2) shows a line schematic diagram of a single degree of 
freedor. integrating gyro. As indicated in the figure, the main ele~ents 
of the device are the svoinning gyro wheel mounted in a gimbal, a vis- 
cous damper which acts between the gimbal and the case, and a signal 
generator which generates a signal proportional to the deflection of 
the gimbal with respect to the case. In operation an angular velocity 
of the case about the input axis (Fig. 2-2) causes a torque to be 
applied to the gimbal bearings which acts at right angles % the gyro 
spin vector. This torque will cause precession of the s>inning wheel 
and its gimbal about the unit output axis. The viscous damper acting 
between the gimbal and case opposes this moticn. The resulting angular 


velocity of the gimbal with respect to the case is proportional 
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boone onpulay Velocity of the cewe with res:ect to inertial-space 
about the input axis. Hence in a given tin2 intervel the gimbal wild 
be displaced with respect to the case hy an angie which is , rorortionel 
to the displacemert of the case with respect to ine tial space in the 
Sam interval. The sional gener’tor sroducss a signal prororticaal 
to gimbal displacement ani therefore its output is pro:ortionai to 
the ansular dis; laceanernt of the grro urit with r=s:ect to inertial 
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POPUL wes Seer Velocities Carries ths an ule) Toren. im Sector cf 
the s: inning wheei with it. This produces a carronert of angular 
momentum terpendiculer to the spin reference axi- of tne zyro unit. 
An angul2r velocity of the case about the spin reference axis will 
then also cause precessier of the gyroscope wheel. The output sizral 
Of tie 7p, ro mal) therefore be in error ty van amount enoport: ona. vo 
the product of the sine of the gimbal angle and the anzular velocity 
SoOuL Unewspin relerer.cs 2x15. 

That is: 
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where in (2-6) the error is expressed as a deviation in outvut 
Signal. The constant k} ie 7a seconscant. of proportionality and the 


angular velocity of the case with respect to inertial space is, as 


Sndicates by the subscrint. aoout the goin reference axis of the erro. 
J 9 s 


unit. 
In addition to the signal deviation producel by the angular 
velocity about the spin re.2rence axis a further deviation is intro- 


duced due to angular acceleration of the cace with respect to inertia 


space about the output axis. Ar:euler motion of the case is transmitted 


to the gimbal and gyro assembly through the fluid of the viscous 
damper. The viscous damper ap lies a torcue toe the zinbal which 
Reseats Mouton Of the-cim sal) Vltn rasp ec, to tie cas2.. Li The gyre 
wheel is assumed non-spinning so that an: torque due to it does not 


acu.) tne Torque “sOnmaticn cn the -ginoal becares: 


@ 
=eC A 


ey : d gir. 


: -~X 

gim ~ gim (I-ca)OA 

where I i is the morent of inertia of the gimbel and cyro wheel 
ee 
about the output axis, A is the enrular accelemation of the 
(I-ca)OA ~ 

Case with respect to inertial scace aocouy the out t (axis sand cy is 
the damping coefficient of the fluid in the viscous damrer. (2-7) 
Simply expresses the moticn of the gimbal assenbly with respect to 
inertial space under the action of the sinwle toraue oroduced by the 


viscous damper. It may be rearranged to give the first order erauatio 


in gimbal angular velocity: 
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From equiticn (2-8) it can be seen thet the effect of a uniform 


anvular acceleration of the cas2 with respect to inertial space is 


to oroduce, in steady state, a constant erngulor v2iotity cof the gimbal 
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with respect to the case. (The transient reriod of this motion can 
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angular velocity ie taken 23 unity “non (2210) ey Oe ae ten: 
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ihesexpressions (2<5) and (2-22) are now identical in form. 
They can be combined to give the total error in body orientation as 


it is represented by the gimbal angle. However, before expressing 
this deviation a word of caution is in order concerning the quantity 
Ream! In succeeding discussion the reference anzle for the system 
will be taken as the gyro gimbal angle, (4,.,). 1s was mentioned 


above Ji stine ro sena) fivit froc Gieith 221s aa ar voloci tr tc 
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er Cwistences Une total deviation 22 tne body 208s trem te inet ual 
reference fram after the eZieccs cf aynaric 12-5 have been removed 


becomes tlc “Combination of (2=5) ane (2-22) 
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Im’ the Gretedine varacrisis t6o“solrces Go error 2 .i1cn ican devo, 
in a system such as the one shown in Fis. 2-1 hav2 been discussed in 
Some detail. The errors have teen compared, and a total system error 
due to bot, effects has tein exzressed. It has teen rointed out that 


the equations of motion describing such 4 system are non-linear and 


hence any attem.t to determine the magnitude of these errors sust 
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resorteto mumeric- le retno 1s or Solution... THiS Tone thal a sascits 
evctem must be cuasen. [he choice of such a system ani th> solution 
Oletie ecuaeOnsr oO! TOvl0m TOrdinis ester. Co ine- aie lo computes 


to det-rmmins the mageatude of the errors =xyre°sed iz (?-7), (2-12), 
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1 = Distance from missile center of gravity to motor 


pivot point. 


Spes(e yA 


ae =e Integrating gyro signel, to motor gimtal 


igu 
angle sensitivity of the pvitch control system. 


S ; ta 
pes(e ) = Pate gyro signal, to motor gimtel angle 


regu? ing 
Serigitivity of tne pitch control sister. 

Apoendix B presents the details of the mathematics leading to Zouation 3-3 

for both pitch andi ay control. 

Gyro Package 

The gero cackase contains three Single degree ci freedom intesracins 
ovros and three sinsle degree of freedom rate gyres. Fis. 3-2 snows 
how these gyros are ::ounted with resvect to body axes. Input, srin 
reference and out: st axes are shown for three gyros mounted on a bese 
which is fixed to body axes. For the purvoses of this discussion the 
orientations shown re;-resent those of both the rate and irnteer-ting 
gyros. 

The rate gyros are assured to b2 cronortional mechanisms. Angular 
velocity of the missile about the invut axis croduces a crosorticnal 
voltage from the gyro at all times. The sigral from thes eyres is 
then of the form: 

(3~4) 


“regu, ; 5 gute) Be 


Where: “rgu(W,e) = The input axis velocity to outcut voltage 


sensitivity of the rete gyro unit. 
The integrating gyros are assumed to be free of drift effects. 
Otherwise, the units are assumed to be identical to the unit described 


in Derivation Sunmery 2 of reference (2). Appendix B derives the 
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verformance equ: tion for the interrcting ryro unit in terms of the 
pyro axes given in Fig. 2-3. Thes2 are then chanced to pive the 
Sauations im terms-of btody-axes com:onents for th: orientation of 


the gyros shown in Fig. 3-3. The result for the X-boly axis (roll) 


becomes: 
(3-5 
ee e ( e 
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“selh eins?) S1Gero Sis Coro get al “nels to %oltome 
sensitivity. 
Sei eo CAs are. ven om the Olnes lye. 6. a Tr eo a 

The eauations (27-1) throurh (3-"), zlone with similer esrressions 
BO Cove soa) “body 296s oor fete. ie ite Une Pet torancte 02 che s7s vem. 
Thes2 might be cor.bire? ir 3 straliechtforwerd manner to zive a single 
performance equation for each axis of the :issile. VHowe2ver in 
prozramuing for the digital computer it is much sim: lor to treat 
the comvonent performance equations as given esove. Tnerefore no 
overall performence equation for the ax2s ar written. 

Appendix 5S sives a collection of the terfemmance ecusticns for 
all components for all three body exas. It is conveniznt to make 
Ceylain changes Ol. nes Variaoles 15 )uhesessonateens at. this ooinu. 

The equations of motion for the missile (3-1) about each axis 
are divided through by the moments of incrtia of the missile about 
their respective axes. This changes the units of the equations from 


those of torque to those of angular acceleration. The ratios of restoring 


tcrjue and interfering torque to the moment of inertia may then be 
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defined as new variables representing the engular acceleration of 
the body due to these toruues acting alone. 

The definition of ansular acceleration due to restoring torque 
can row be carried into the equations for tne torque generating syster 
by dividing these equations through by the apprenriate momants of 
ineGtia. <inils ~ubs both the body and torcus pererating system Scualions 
in terms of angles and their derivatives. 
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Where: Scie ) = Gyro gimbal angle to restoring angular ass-ier- 
gim’ m 
atLoOnesens| Uiviey of the Soll controls, steu: 
5 ure) = Angular velocity to restoring angulzr acceler— 
res(W?"m : 


ation Sensitivity sof ta] roli e@antrol system. 

It will be noted thet the equetions (2-2) have been modified by dronpins 
the subscripts indicating body ax2s compone:.ls. This is in agreement 
with the system equations of motion, Angles and their jlerivatives 
wiich are referr2d to the inertial axes set are identified by the 
subscriot (I). The sensitivities defined for the roll control system 
above have identical ccuntar parts in the pitch an? yaw control systens. 

The constants required for the set of equations (3-5) must bs 
sp2cifiea.. The comblete 3: of these constants for the -sys-ex are 
given in Table ITI-1. The inertia ratios were arbitrarily chosen. 
They represenl what ‘aight b2 found in a missile of small *o mediun 
$1ze such as the second stage of a two or three stage vehicle. The 
integrating gyro characteristic times =r2 thos» for the MIT Instrumr.ta- 


; 4. : P 
tion Laboratory, 19° -crro. The natucil frecuency and damping ratio for 


—_ 


the torque generating systems were arbitrarily chosen. The values 
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chosen are representative of what might be expected from a hydraulic 
servo in combination with the engine inertia. The system sensitivities 
which have al] been combined in the torque scnerating systems remained 
to be chosen. These were chosen on the basis of é lirearized stability 
analysis for the system which is discussed belov. 

In an effort to obtain sore infornation on system stability, the 
equations of moticn for one axis were combired., All coupling terms 
from the remzir:ty.g axes were assumed to be zero. The root locus of 
the characteristic equ-ticn of the resulting fifth order system is 
shown in Fig. 3-3. The calculations leading to this locus are zviven 
in Appendix 38. 


The locus of Fiz. 3-2 shows two oscillatory modes. As the oren 
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loop sensitivit 


the torque g2neratirneg system moves toward the rieht half-rlane. At 
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the same time the pole vair et tne origin gue to th Cerraving © 
and the missile performance 2quatior move to the left. ‘he system 


open loop sensitivity (5S ) was chosen so as to give both 


pes(As sd) 
pole pairs the same closed loop dampinz:. 

The preceding paragrephs have served to develop the perfomnance 
equations of a simplified vorsicn of a geometrical stabilization system. 
It has been pointed out previously that the errors under investivation 
occur regardless of system dynamics. This fact can be demonstrated by 
varying the dynamics of the system just developed and solving the 
problem for the seme interferirg input. Since the problem was solved 
on the digital computer the variation in dynamics can be acco:.plished 
by rather elzerentary chanses in the computer prozram. Therefore another 
system is proposed. The body equations of the s3t (3-5) vere retained 


but the torque generating systems and intesrating zyros were assumed 
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to have no dynamics. A set of equations similar to the set (3-6) is 
given for this system in Appendix BP. The programing of these two 


systems for the digital computer is taken up in the next chapter. 
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CHAPTSR IV 


METECD OF SOLUTION 


Non Linear differential equations are not in general amenable 
to soluticn in closed pon Finding the exact solution of a number 
of simultaneous, non linear differential equations is ever, less probable. 
However, numerical techniques have been developed for obtaining partic- 
Mlarevsolutions of diiferential equ-ticns, cither linesr or nen lincar. 
The soluticn so obtained will, as the name sugvests, be in number form. 
FOr solution, the equations must be stated with nurerical coefficients. 
The coefficients may be variable, but the manner in which they vary 
must be known. Certain initial conditions of the variables are also 
needed. The computation time recuired to solve for the dependent 
variables over a given range of the irdependent varieble, will depend 
upon the numerical technique and the accuracy required. For imvroved 
accuracy the increment of the independent variable betweer calculated 
points must, in general be decreased, thus increasing computation time. 

The :.umerical method used in solving the set of ncn linear 
differential equaticns (3-6) which are listed in the previous chapter 
is the Runge-Kutta Method with the Gill modi tucatien. Use of this 
method requires a known value of each variable and the first derivative 
of each, at the start of each integration step. If the differential 


equation is higher than first order the integration process must be 
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used as many times as the order of the equstion. For example if 
2 


oy = a ra N Grn Ae @ espe ) the value of the second derivative 
dt at eat 


of y is calculated at a point, say to» where values of x, y, théir 
first derivatives, z, and t are known, The integration method will 
give the value of the first derivative of y for t = t. - At where 

At is the integration interval. To get th value of y at the point 

Pe + At the first derivetive, whose value was calculated for t = t. 
from the previous time step which started at t - At, is integrated 
using the same process. The equations wnich indicate the need for 
solving for lower onier derivatives are called aD tary 2euatio 7s. 

The eight second order differential 2eoustions m-thematicelly describing 
the system, therefore require that the integration method be apnlied 
sixteen tires at each time step. The transfer to inertial axes of 
rotation rates excressec in body axes give three more derivatives 

that must ve calculated and integrated for each time stepo. The eouatiorns, 
in computation form, along with their auxili-sry 2quations, are shown 

Im Avpendix Bb. 

Deeds lemuoo! ised sol tet ene Ul wee ren bra bvecua Ul OnSs sol yen 
Mehopendix B Was the 1. 6. bk. 600 map etic drum data—srocessing 
maenies ) 2) Tne unit used is operated 5; the Mathenatics Group in the 
M.I.T. Instrucentation Latoratory. The machine used has an internal 
Storage of two thousand re isters, two tape w.its, floating point 
arithmetic, and index accumulators. 

The problem was programred using a floating address ass2mbly 
Routine. ‘called Flaa‘15) developed by the M.I.T. Instrumentation Lab 


Math Group. The basic advantages in using this assembly routine are 


the simplifications in coding and the optmizaticn of the program to 
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reduce computation time. The differentia] eauation routine was 
available as a subroutine requiring only entry and exit programming. 
The program could te broren do.m into functional parts as follows: 


1. Obtain the forcing function, Ae for each axis. 


int)? 
2. Compute tne values of the derivatives anit store them in 
the proper r2gisters tc 3rcuer the 
3. Sntry end exit instructi:ns ne2ded with the subroutine. 
1, Put tre-resuine 7 mame cs 1 the. ces. rec form rants read 
tiem On Of the c: mpnter. 
m. otart tne. cycle Over esata et o> sex. timesster. 
imeranhiin= the vrodles te cite se=r nezded to five the ssouired 


L Sas 


accuracy of solutzon prove: to be unfortunately cma. d ihe root locus 


Oiwune Mineerized system, .21Von ih Une Orece sting Charmer indicaved 
eet ties DeSuulat-d Sysocm Ssnoul. ce 62acls./)in verging cite the 
OrOsran 2) Sinusoie! ©orsin= funetic:; of 15 eadmmes oor seccnl 22s used. 
in a linecr system time increwents of one twentielh of a period usually 
give accurat> results. For such a time ster the sclution of the system 
2guations diverged. Using one millisecond time steps si:ow2d that the 
solution tas stabie end so=ue of the varissles charge: rapidiy recuiring 
that. the time ste = te smelt for the dicital solutior te: follow the 
Cragin “J55 (6 yg enero riot at Ge yoo sea eo a LPeoste C1 ave 


millisseonds the com outed values of the varinbies agre:2 with the same 


12 seccr i was used. “Atrten -ciliisecona 


) 


comvuted values «wien one ni 
time (steps Cie) solution dee ersee. ihe rercre dive mitivsecond, tif 

steps were chosen. This time step resulted in approximatly twenty one 
hours of corputing time bein= required tc simulate one mir.ute of problem 


time, using squar2 pulse forcing functions, 
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Suggested methods of forcing the equations werz by the use of 
actual acceleration data, summing sine waves, anid using pulses. Use 
of data would require reading data cards at each computation point. 
This would have recuired the punching ofat ‘east four nmmired cards 
for each seconi of problem time. Evaluating a sine takes 150 milli- 
seconds. For a single sine wave per axis almcst a h-lf a seconi would 
teretcre 52> recuired each tims toe forcing function wes calculated. 
Porea squars pulse onl additions and subtractions are rec iired ‘arid 
lice Lie Toss ont 1 he orcine funciicn bart ol tne procram. ihe 
Beleevlon Of SQU=Te TUlSses Uses wr this Anvset oot sss oe Primarily 


determined by the conpul+tion tine reouired. 


Ther syster was recramncd. CO. a2UP i. ths reanons- tortnre 


ie CAM) Wi) ts) Oise ectrent Of the £1Ssi2e2 abou. 31 three 


aAXCS. 


3. Random rectanzulir puls?s of angpilar ace2zleration acting 


on all three axes. 


ssile atcut all three ax2s was 


pe 


the r: 


be 


Whe [aitial  aisplacswmenit. 6 
us2:d to determine the transient vehavior of the system. A similar 
iyo Amish -Ciscoiee ent cCondiuion cold o> Frealiveld “in an actual 
misstle. For a multist222 missile the sey aration orocess would be 
exnected to ceuse some rotation of the next stage. ‘when the propulsion 
system fired the control system woull see the disniacement as an initial 
error in orientatior., 


Rectangular pulses were used for 'eriodic forcing of the system 


(46) 





primarily due to the speed with which they could be generated, as 
noted earlier, Use of these pulses added no appreciable comrutation 
time over that required for the initial Jisplacemer.t. case, These 
pulses had a maynitud> of four tenths radian per secon! square. 
For each axis the pulse duration was a tenth of a seconi with two 
tenths of a secon! between the en? of er :uls>* ant Lhe stest of 
anyviyes, The pulses were phased so only one axis at a time was being 
forced,, “wher theipulse vos taker or: of one axis and aprlisd to 
aaCviervits Sig was Cadmeet, iis sev ehe ;erioa of the lore ne 
function's fundarental frecuenc:’ as six terths of a seccund. 
The random culses of angular acceserati-n -ere limited in 
Happier aero tO nine mol lit iiens cer second S71 rod in 
lrcr2ments of ten. ruls2 duraticn was limited to five hundred nilli- 
SeCOnG opal. Th Crome) bs 6. five fi lisecerniS.6  hierslen. man. ude ene 
Cuneo Or Van OUISe Was Co Lerma Sse Usin = Gie1ts i <noreris ters 
containin= calculated valu:s of scnm2 vari:sble. These rezisters 


GCemvained Valucs in) 22¢atin= -o:nt Porm. ocdizits loss tier txo 


pee 


digits to the right of the first s .r.ificant figure w-re used. These 
$US Varied at every time step... he digits used to dever7ine the 
Porern2 Puneuion os cee axis can fren Varia,2¢s orm rds ini enced 
by metions about the other axes. A singl2 digit determined the coulse 
Mapiavude, © .0 determine tie Sign Oo: «Ge DUls=, five MMs subtracts) 

from a Sinsle cizit and the sign of the difilerence was usei, The lereth 
of the pulse ir. five millisecond incremerts was determined using two 
digits. These procedures give a forcing function on 2acn axis which 


is independent of the mode of forcing on the other two axes. The 


couplirg between axes given by Suler's ecuaticns ard the gyro equations 
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does not however maze the angular motion sbout any cne axis indenrendent 
of the motions anout the other two axes. 

Progranming for the further simplified systems which were mentioned 
in Chapter III was accomplished by removins parts of the original 


program. 
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CHAFT“R V 


DISCUSSECh Cr Robi 


This charter discusses the results.obtaired from the solution of 
tne system equations on the digitel] computer. the discussion 1s concerned 
Primarily «in tec ansilar ditferoress cevelorpe! in Chegusr Jl.) liese 
Git {erences feureseti. Srrors, lm imei cation 0: body oriertation. he 
difference <2tween icertial angie and sody anvle is the error in 
Orientaticn due te finite retations. The difference b tween body angle 
Wie eyrO pli Dalsaneie 1S tyecerror in Orcen ib «nm cause: bY 27710 
characteristics. The sum of these differences (A,-A eis. the en 
errer in Orienteticn due to these eifects. These differences do rot 
eneiude Che -asVielion an orientation resutvine 2roniimaeic lags in 
Rie S55 Cem. 


Dransioaietes nonce 


. 


It was desireable that sone examination of systen stelility be med2. 
inorder to do. this the resvonse al the system tov~an init .a) ancular 
displacement was computed. io interfering angular acceizraticns were 

n.Jied durin, this response. 

Figs, S-1, 3-2, and 5-3 show the response of the systen in restoring 
from an initial angular disrlacement. In these fizures the gimbal angle 
is taken as the reference variable. This was done because the gimbal 


angle is the variable which the system attempts to nul]. The figures 
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for the pitch and yaw axes show a demped oscillatory respcr.se. The 
response in roll was only run until the gimbal angles on the other 
two axes were essentially zero. St-bility of the roll axis system 
was not a problem sirce prorortionel control was used in th: torque 
generator and the loop sensitivity was relatively low. 

The root lcecus of the linearized pitch control system was discussed 
ipecneapies Lil and »showno2) Picy 3-3, The close loon gait, “pave a 
damping ratio of 0.445 for the two oscillatory modes of the system. 
The control system in pitch and yaw is a fifth order system in angular 
Velocity as siovwm. in Che te~ Tile ihe systemis also non linear, cub 
some feeling for system characteristics may be obtaired by comparing 
it with a linear secend ornler system. Charter 19 of Ret. 6 cortains 
curves of damping r' tio versus transient peak ratio for linear second 
Ord2r systems. . 2ntering tnece curves with tae trenslent teak ratic 


obtained fren the system response fiv2s 


RY 
{i 
¢\) 
:3 


Nping ratic. a denmred 
bere Ol Osc vllaticn may aiso be obviined from coe traneie. resi 0:c2. 
An undamnped naturai teriod may be calculated using the dempine ratio. 

mae ean eqene Pes ous]es SnGwn “in Pigs, S-)° and S=2-are nev those 
of a seconi Jrd?r s3stem. Therefcre values of damring ratic and 
undamped natural frecuercy b-sed on second orier res onse characteristics 
are m2zaningful only as a rough measure of the system characteristics. 
From Fiz. 1, the yaw axis response indicat2s a damping ratio of 0.48 
and an undamped natural frequency of 1.19 cycles rer second. From the 
Pitch> axis response, fig. S=2%a dancing ratio cf 0,461and-anm undenesd 


a 


netura) frequency of 0.95 cycles per second were obtained. In Fig. = 


ee 
a 


the slight recurvsuture of the gimodal anvle curve indicstes that the 


roll control system is rather heavily damped. The transient response 
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indicates that the system is dynamically stable on all axes, 

The differences between inertial, body, and ginbal angles are 
plotted in Figs. 5-1, 5-2, and 5-3. No interfering angular accelerations 
were apvlied to the system during the transient response. The existence 
of these differences during the transient response indicetes that any 
angular mction of the missile ma: cause errors in orienteticn. The 
absence of any means in the system of sensing these differences is 
shown by their constent value once the gyro gimbal angles have reached 
their steady state nul] bosition, 

Information on the source of the differer.c=s can be obtained by 
consideriny the orien'ation of the gyrc axes relativ2 to the missile. 

To aid in comparison, the body axes orientetion of the gyros are 


tabulated below. 








pS Gyro Aes 
Rotation Sensed, TA SRA 
! Roll X Y Z | 
Pitch 3 x ~-Z ! 
Yaw Z Y ~-X | 

Table ol 


Wee Ses LLG Groot AT ICN 

In the equations (2-13.09) fo the integrating gyro, the terms 
introducing differences between body and gimbal angle are angular 
accelereticns of the missile about the gyro outvut axis and the product 
of gyro gimbal angle and angular velocity of the missile about the 
gyro spin reference axis. 

Comparing the body to gimbal angle difference curves for the roll 
and yaw axes shows a marked similarity in the oscillatory character of 


the differences. Both gyros are oriented with their spin reference 
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axes along the pitch (Y) aais of the missile. This would indicate 
that the dominant cause of this difference was the result of angular 
velocity about the spin reference axis. The greater magnitude of the 
error in the roll axis plot can be accounted for by the fact that the 
gimbal ansle remains lerz2r for the roll axis than for the vaw axis. 
The relativel, small effect on the body to cimtal anecle difference 
Pesultimesirog acpular accelerations acoul tie ouvcuy axis is indicated 
in the vitch axis plot. The :-itch axis ‘"yro hes its srin reference 
axas alone the roll axis, . ine siow rate of chars2 ~of 2izo4] anele for 
the roll axis eyro shows that the roll rate of the missiie is low. 
This indicates that the :iten gyro is little affected by angular 
Movavi Or ALOUL 10S Slit, seis rence axis, §ine ouvcut axis sof the: o1ten 
Pyro 15 along tne yan 04) exis. MW ohe-curveture co Se Z—-2xa 6 S770 
ginbal ansle curve i:.dicates the nagnitude of the ar 
6: the missile about tnat axis, It will be 2:oted that the peaks of 
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In this transient resj;onse the ansuiar acceleretion of the nod: ebout 
aueyoo OuLpuL 2knsS ageears to resule a Jictic error Henever ne 
general statement s22mns warranted for the case of arbitrary missile 
motion, 

Tne inertiai to een arcle dit erence "curves ers =lost inentical 
for the pitch and yaw axes, while this difference is almost nene>risternt 
in the told Axis DlOt.) casas Gaflerence asrths Tresulu of iviterota— 
tions as was shown in Chapter II. ror the roll axis the difference 
2quation is: 


(5-1) 





Bore MIC Ne Canes ene drs rence sS likewise the deffercnes betwen 
the products of .ody anrle ani angular velecity comronents associated 
Withiuas OU,cr UNO axes. 

Comparison of Figs. %-l and 2 shows that the gimba] angle curves 
enevaltostsideitical. “This means that at any time during the resronse 
tne LOcy anegle: ani rotation rites ate nearly identical for the Y and 
ZeoxeSe. The tivo terns i. the difference 2zquatior above therefore sun 
tO aamost zero. “nis also mares the pitch ard yaw axes differences 


a e e . e a e 
Basco ri@omin cel wathotis* 95 Goon. ir the Pfivures, 


Gieprstaye Of 2 tlw le Sstae> sissiie.. tn the 
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thse Ml tial CON etl Given Une Sol ero) gosto a. tne, thristren.1 se 
Pore sens “CONnLrol. 1s commenced, = If uUmevancule> velcc2 tiss ore smeaid 
SndeCar i. 1Lecescted, this wouid 2stecdish 2ritial conditions similar 
to those used for the transient res: onse. The body uwculd have a 
GSveat2On in Ofer Calicr Wi the cero solztin. ratesmy Tor the, response 
Salecuiated Voevit tere.ce  e2tween icer ie in tie ern body angle wes 
eo roxiratet) cre ver cent of ths inttial distlacemert. for all thse 
anes f 19s diiiscren] Os the tots] sonientst.cn error deve lop-o (dumine 
the respor.se. 

fae Inlor att nN Obtain. on Tals calculate. ol -viee erans eit 
respor.se of the system is as follows: 


l. The control system is dynanically stable. 


2. Rotations of the missile rroduce errors in orientation 
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due’ to the eitects of three dimensional finite rot-ticns -ni th: 
characteristics of the intecrating fyro. 

3. These errors in crient:tion are net sensed by the 
system and exist in sterzdy state. 

4. For the response compute? the total steady state error 
in orientaticn was aivroximetely con® «er cert of the initial displace- 
Mets. 
hemiodic Porcine 

The systens resverse towa veriociec forcine function in interfering 
aneular acceleration was obtained. For reasons >lready discussed in 
Cnecter IV all of the forcigeeitici ions use? in this analysis were 
Goma sede (CC (ONG Oo Segue] Culecs,. PCrvcn= ser iodice forcing 
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Pimclicn: tie: S22 o) Of ls > s.siowniim Pies) -4) eo and 5-6 was 
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Oreicraridy.sej2ctaq “One “22 Vetieex1S JO! (Ure ssien eee pulse. 

ére elternatety positive and nesative with a period of six hundred 
Prlilissconds.)  2acn ose “as a Sureti cre! cre ii: ec i issconds 
and a magnitude of four hunired rilliradiers ner seconi scuared. The 
ig of the pulses cr tne various axes wes arbitrarily cnoser: so 
Cheat che Yaw axis ses tc Uo irs (Ol sa period bening tne Titew exis 
and the roll @*is was two (tnirdseol a pericd benince the ‘yer axis. 

The veriodic forcing function just describes 1s equivalent to a 
linear Vibretional acceleraticn cl tne missiie~of -four Tse. cer -s2ccrid 
squared at a distarce of ten fest from the cer.ter cf gravity. In the 
abser.ce of system restoring torque this acceleretion acting for the 
duration of the pulse would produce a displacement cf about a querter 
of an inch at ten feet. 

The response of the system to the ;eriodic forcire function jus 


ve: 
re 
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described is plotted in Figs. =< through *-9 for all axes. These 
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piots show the forcing functicr, the gyro gimbal any]? end the three 
eueviar differerce= for sae, axis for the first ten seccras of system 
operation. At tne end of ter: secunds the system Lad reached steady 
State on all axes so that further conputation was unnecessary, 

The ;lot of the syre gimbal enzle is essentially the same for 
all axes. After an initial transient pericd the gimbal anvle settied 
down to a steady oscillation about its rull rositicn which ancears to 
be almost sinusoidal. However, cormyparisor of these nlots with those 
Gl igs. =-) Larouzn o-3 shes that. the gir bal anzle res sense is 
actuzlly maje uy of a continuous series of expcnen tial curves similer 
to those oktaine: for the traensicrt respor:.te,. 

The behavior of tne gimcal angle during the initia] transiert 
ceriod is consister.t with the results of the trensiert anelysis 
discussed :revicusly. Zach exis uevelops sore dynamic error during 
tne transient period uhichn 1s ever.tusily damped cut. In the cas: of 
the pitch and yaw axes this darn:ing recuires arproxinately tic secondcs 
while for the roll axis approximately six seconis are rewired. 

The effects of courling ermcnz axes in tne systen. can te seer: ir 
the riots of ginbal engle for the pitcn and roll axes. sSacn of these 
axes 15 initially unforced and “et in =ach case some small gimbal 
angle develops immediately due to the motion of the missile about the 
yaw axis. 

The angular differer.ce:. of Fig. S-4 through 3-9 all develop in 
a fairly consistent pattern in spite of the verietion in dynamics 
between axes of the system. Some general discussion of the basic 
patterns is warrented before considering the results axis by axis. 


The difference between body angle and gyro gimbal angle is, in 


each case, an oscillatory curve which slowly drifts away from zero. 
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In steady state this is a combiratior. of a stesdy oscillation about 
zero and a constant drift rate. Tnis can be expleined in terms of 
the two error producing effects. discussed in connection: With the 
integratin- gyro. 

In steady state the components cf anguler ecceleration of the 
missileare symmetric about wero in the present. cese. Therefore for 
@ full cycle the net ar.gtiar acceleretio:r cf the nissiie ebeut e 
@iven axis is zerc. The difference betwee: body and gimbal angle may 
Scevrlate cue TO (the elfsct ol anvuler %cecle sation, abveut.4 fre 
Gucmuteaxis but for a full cycle the net dil l@rerce creduced Uy its 
etvecriis zerc, <Any Uriftoim the :dot- af cay te sintal an<le 
Gisterence tustethen te cue to the effect of ansuler veiocity of the 
missile about the gyro sj;.in refe.sence axis. 

The product of gyro gimbal] angie and missile anguler velecity 
about the gyro siin reference axis wili not in generel irteg:ate to 
zero over a full cycic. Therefore over a full cycle a net difference 
between body; and gyro zgimbal engles wil} usually exist. In steady 
State the diflvererce due to this effect will be the same for each 


S the linear increase upon which the oscillator; 


Lb 


eveie. [his predu: 
SlieCUsS ere si riper . 

The carves for the difference between inertial and body angles 
ere characterized by a small oscillation, the mean of which is lineer 
With w.tanite slope iy stoa3y Gtate. This linear cavervance 15 to ike 
expected for the case of periodic forcing. A net difference deveioned 
during a given cyciz of the forcing function would also ha develored 
in any other cycle in steady state. Therefore there would be a linear 


accuwiulation of these differences with tine such as is shown in the 


Curves. 
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The sum of the differences just discussed is the 7ifference 
between th? inertial ar.22 and the gyro gimbal angle. This difference 
is also plotted for each axis. It represents the total error between 
the reference orientation and tke orientition to which the system is 
Currently iattemptigie te drive tagimissi ie tue. to all-of- the effects 
under investivation. 

The angilar differe:.c:s for each nazis at tiie end of ten seccnds 
are given in table =2. Since all Jictfersunces slerdiiy 3:.c-ease in 
baenitude: these ere the aiaxitum °roors inert oion wich the 7s stem 
experiences for the ten secon) run. ror hOS2. 11. Miers: c¢>s +hicw-are 
Soci livilory tne mumber>s Piven ere cle sea error Ac ter. seconds, 

Lae. ore S-Di Mee On tie Sale os 2 Pe peice cee ae 


mccain SOL LUNG “Or Mie Sst less iar Oosci dd Cie Wie stend (oe -to, 


* 


= am 7 


Dimes Gis vrO SSLSi Livi, Irom GOi 7 anei= Ceo wyre “17,041. =" on) 1s 


3 4. oe lc | *~ . 73 ‘ wy 7+ ~ vA a > =) . R aie + eK 
Ceo, se ee SCAM? O. SyTO il OA ld temrte Ose tation cen ae 
Ueed. 2m lOmnin= set1Os 9) Oe Very Le 2S rere T2asin Uk Son redets 


+e. 


SLrOrs JtNcOFi tics clon Lo OCs, a lor 1.0L lon. “lu nsnok sa) ce ae led yas 
the magnituce of the body angie is 1.t tunat of the gizrbal an-.12 due 

bo the errors, “HowSsver. (Une ei ude Of Une Dod an 1S OSC iia st ions 
are the same as the amvulitude of tne si.esl] angle ose ilaticns. 

Tier Cees ceri Of IFimininS “Gime sl(Or wiles vie suits Of Senlodic 
forcing have been i2ztermined i¢ sncert compared to the firing tine of 
most missiles. A good approximaticn to missile firinzs tim2 would :e 
ore hundred seccnds which would reuvuire a comruter rur. ten times as 
long as the cne that was mad>. However, once the system has reached 
steady state operation and remaired there lons enouzh to establish 


the trend of the errors further com: ut~tion is unnecessary. The error 
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CURVES seccins linear stay bose <tra pl lst to any lssimpe value.) In 

the present case this extrarolaetion has ben carried out ani the 

Re SULuones errors sm Orient to. ate snown in calee Te elorng “ith 

the values for tan seconds. The vaiues of the differences at one 

nurdea “seculs are (sived ir Voth radians and dep-ces for converierc:. 
The follovwir.e informsati.r. «as ot tained from the resnonse of the 

Sveten UO a teri Olic- ittterioerine arrular 2 -celeratvicon;: 
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Cores daronsiona! fini 


Oe LheSe ce ie Gre a Seventeen im “svecd: State ited 
Bay es VOU wer se I won ence 1 Ole ee yealGe orto or 
Hes COnes tas 0. Jere as ity Sole OG mi) ave aS ep bon. 


ORO oe? ad haed atone is Buea dere S el ore ee ae 


MGeor Une an diwude er oSsine Snela> wocien-about, Cle Sxis in-cusst.c.. 


So. IpiliGe Une ie rease DMN eerrel ts Janeen se iraeco lac tom US 

190 seconds was .ossible. 
Cie Ust, 
Random Forcing 

In ord]r to simulate imterferin 
realistically as vossible with'n the limitations already imposed by 
computing time, the res:onse of the system te a series of random 
puls?s was determsined. Since no steedy state could be exrected for 
such a fcerceing functicn the probl?m was run for forty seconds of 
real time. This required approximately ten hours of computing time. 


The random forcing function was rade random «ith resnect to 
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Libcemotybie pilsce cn. tacts rstleGs. ine (Sipnele Gere ules = vs 
Mopac. dil. te duration were ind?pend3ntiy caosen ine random 
fashion. This was done automatically in the computer. The sign of 
the pulse sas chosen first. Th mechamism:; lich this was done was 
tive slpas Sch Gal A robsbir lity or 


7 


ar) 


Suc that OOSlCiIVve and 12 
eccurrence Of one half in a civen selectian. The magnitude of the 
SUulses wactlinites to the values fron O to 9O mi liiracians ver second 


Of Ueivan: livras ane e> secon Ssausered.. gach 


UY} 


Ssuireis in stec 
Value had a provabiiity of occurrence of cne tertn in a ziven select- 
Pon, lne -tuivat@cn.or tha pulses was se. 20 24 *:t. ramen, in the sae ee 
EvoOnrmo-toeoO) mil. Lseconds Sith selecti.n- bane leit 
SWerivermtelisscOnis, ine -Torcite 2unccicn: Clcgitie In Chis manger 

es ISO WMaOr Sacra tis ty TleS, 0.1 Oe roach ol aw et ould. be oe 
VeriOuseaxes Since sélection on each axis was como letiy 24> endienv 
Omtie Other twee. 

The response of the system to the random forcing function Just 
deserieced Ps SaO.msin Fizs.2 10 throuch) 7-23, 

The system response as siown by gyro simbal ansle is about what 
could be expected in view of the previous results. The ginbal angle 
oscillates randomly in essentially the same marner as the forcing 
function. There is no steady state condition arp parent in the response, 
None was expected since the forcing function is random. In general 
the behavior of the cimbal angle is consistent with the rrevious 
observations on system dynamic stability. 


The angular differences have the same general character as those 


observed in the periodic case. The body to gimbal angle difference 
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CUrvesis ‘oscillatory, while thoe-inertial to boldly anzle -lifference 
is relatively smooth. However, there is no systematic divergence 


Ol these anmilar ditferencss Con tnis run. In this resrmect. the 





response was different fron the results for veriodic forcing. 

Since the response to random forcing is essentially statistical 
in nature generalized statements cen not be made on the basis of a 
Single run. The anzsular difference curves shown are just one sat of 
Slipime lite Mueder “iat Could Occur as a “result of This tise of 
Porcine. [he fact that tne ansuiar ditierences Obtained du-ine ons 
Pltingo COU Vero: @acrs Holl trolly thateta-y would Act diverce in 
ame ner run. woVever, tie Gi fle run snowed es demonstr le Unat 
Grrers 1m Crientation db exist Sver wach. c2riom forcirz... “che 
hatereerine an ular ecce leet On am ULS noes -Ol- ce ferioiic for 
Wiese Gricrvstionel 2rress tocdey2loc3 und “Un2se Are oie lords conc— 


PUSMenGebl encalt TaoOne Ol arr un dion Liss CS cole Omer ama 


Gneraliaed ressits WoOuls St. be or aiid. 

A run was made to show that orientational 2rrors will exist 
regardless of control syste; dynamics. The system dinzmics were 
eitered by assumin= oropertional rescons> inthe interratin=«: evres 
and the torque generating systems. The reduced euuations are shown 
in Appendix B. 

This Simplified system was forced by the same periodic forcing 
function used previously. The response to this forcing is shown in 
Figs. 5-14 through 5-15. 

It will be noted that the plots of this response show missils 


body angle as the reference variable. For this simplified system 
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the gyro angle is identical to the body ancl >, Therefore the only 
difference that can occur is that between inertial and body angles. 
The angular difference curves fron inertial to body angle irdicate 
the error in orientation develops in the same marner as for the 
complete system. The error is siightl: less in this case since 


thescontro. system rescons: is faster. Therefore the missile will 


(9 


experience smaller engul:r deviations and lower angular velocitie 
Foren Same .OTCing [uCclon. ) 21S a2Cecunts scr Tne redecuion. in 
the net error in orientation duc to finite rotec en 67 facts. o In 
any case the error does exist regariless of control sistem dimanics. 

the results obtained in this investigation srove the existence 
Siecerrors IN Orientaliciias to Unree cine ts1eral fi te rotations 
ANG inleeraving Byro  Chetac  @6iSties,. Some sini Oration “as Octet 
LOr Unie cCages Oly SrIOlle TOvem Sanden ove ae. ato trera cs 
rstomse froin an initial anzgular displacement. 

For the case of perioiie foscine = surther Mives ire Mas 
needed to determine the effects of tne frecuency and magnitude of the 
foreiny Cunctiom. “ln the casé of random forsi1e ore “runs wilt Se 
recuired to rermit a statical evaluation of the errors deve: coed. 
If very many runs are to be mad=, a faster mehtod of ccmputation is 
highly desirable. This can be ascomplished through the use of a 


Pasuer Ccrpuler -sor Usraaps by Charing the Tethcd of Com ita lio, 
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may cevelon in the segue ri cal stabtilic=.1cm =). 
ByerOs.. “Lhese “errors are Jisinct»irom vi t2voreer- in Orientation 
which such a svster would ex.erience Jue tu sistem dynamic lags. It 
was demonstrated tea ee Cm rea (Coe Ons be ie, Src rs sores 
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Those errors G=v-25 2) di] te th= courling of messi l= an=iler motion 
Betvetn 225 (0; “Tae Ss 7sle7. 

when tne Ys,Ster was ferced witts of veriolle interverits nzular 
acceleration o1feil -24S . Unere vas a Systeatic increase in orientation 
errors wit times, “Alter ten seconds tis total Serer itn-orientation in 
yaw wes forty per cent of the suximum am:lituld» of the enzuixzr oscillaticn 
Gtatie Missile ay steady state.) Correspcncing 1-es-icr- tae-oiten 
ane Mollwaixyes were tyenly sever: er cent and thirty? reur per cent 
ReSreceelvely. sakese Crrersvincreasce Jinearly with» tunes. nere cr rerac 
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one humdred sscords the errcer in yaw becomes feour frundired ner cent of 
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tiesto apis ties oc the anegulas Oscillation of the riasile, 
The sam2 kind of errors in orient:-tion developed when the 
PiUerie Pili Se cller e2cceleraulonis sete random. « .en2ver.<iue te the 
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Chic ee Sa CCl lo COLE) ee eS  SarS 


Fie. A=1 shows. the inertiel and body axes sets deseribed in Chapter 
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Ii vane an Sruulrary vec vor Yowhich is fixed in the notly axes set. e 
maou er 6s LOA°*Ursss 2h> eo orents Gf 2hs se “Gir in thevanertial 
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have assum=d an aru_trery orientate. «5 <sier-l, en =rtitrary 
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cf the vector V in inertial ees are derived, tfter three rotstions 
ene results for the @pdivibwl cotaticons are combined to 
COmPOnstivs Of th= vesuGr in inertia: axes in tes ci Une three 

areres if rotetion and the comcontnts in bedy axes. Tne resulting 


& 


transformation ecu: tion is not unicue. Thers is a differert trans- 


formation, fOr 2aceypossieic: Order Ol retebions. “coveves, a2 sSimplatieataen 
US.made at, tne ena which pe-s eround this ditficulty. In the Gerivetion 
below the orier of rotations is about the 4, Y, and Z-axes in thet 
Ons. 

Letethe body aves of fiz. A=-1 be roteted al cut the x inercial 
axis through the angle A, as shown it. Fig. A=2.. From the figure it 
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can be seen that the compomm@nts of the vector V in inertial axes 
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these all reduce to (A-?). The equation (A-9) and (2-2) are identical. 
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The remaind:: of the system performance equations are concerned 
With the determination of the restoring torque (x) as a function of 
missile motion. 

In the missile configuration being considered restoring torques 
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Substitution of (b-12) (B-13) (B-14) into the s;stem equations has 
eivewet fect ol fuiing the syeten sensitivivy for each axis in the 
torque generating syste1. ecuations, except for the signal generator 
sensitivity. If the gyro equations are expressed ir. term:s cf gyro 
Zimbal angle instead of voltaye this sensitivity is also carried 

Mier IMtO. Une oOroue semen ting Syolor, bk Ae ee Cm nee es Cis 


for the systems then beccre: 


® e080 ® 6 
Re ae ALR 1 eA + of 
x ae ay (int), my 
@ e@¢ e@ 
Tf wT <- 
"ay (TR) * = Bea . a 
e ® ©@ eo: 
ee (Ti) ey ~ es ar 
ral ee F n ee i 
a ~ “Cres es £) “ein i“ ees ee) ie 
a pis x 
= z 
oe Sy ta NE Nie 
2 Oe) ae ea ie eae 
7 re: eee mee. a 1) Ans 
RR i ir ofoks 4 @ oo mmeaie yy Me 20 
eee es 4 y 
Les 1S 





=~ e¢ .6°¢@ an 2 
a ~res(3A) “z 
ee . e 
(Cr). + A nus a A = Be — eee a 





rr ae 





5 € *(res)(e 3M) (ren) (ise) 
(res)(i;a) ~ ———-fa4 


x 





and identical definitions for the sensitivities hold for the remaininz 
axes of the system. 

The equations (B-15) are the set that wa solved on the dipital 
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The set of equatiur.s (B-15) was the set solved for the system with 
Simplified dynamics. The constants her: identical with those for 
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